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ABSTRACT 

Prior to the incineration of a white dwarf (WD) that makes a Type la supernova (SN la), the star "simmers" 
for ~ 1000 years in a convecting, carbon burning region. We have found that weak interactions during this 
time increase the neutron excess by an amount that depends on the total quantity of carbon burned prior to the 
explosion. This contribution is in addition to the metallicity (Z) dependent neutronization through the 22 Ne 
abundance (as studied by Timmes, Brown, & Truran). The main consequence is that we expect a "floor" to 
the level of neutronization that dominates over the metallicity contribution when Z/Z Q < 2/3, and it can be 
important for even larger metallicities if substantial energy is lost to neutrinos via the convective Urea process. 
This would mask any correlations between SN la properties and galactic environments at low metallicities. In 
addition, we show that recent observations of the dependences of SNe la on galactic environments make it clear 
that metallicity alone cannot provide for the full observed diversity of events. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: general — white dwarfs 



1. INTRODUCTION 

The use of Type la supernovae (SNe la) as cosmologi- 
cal distance indicators has intensified the need to understand 
white dwarf (WD) explosions. Of particular impor tance is the 
origin of the Phillips relation (Phil lips et al.l [l999). an essen- 
tial luminosity calibrator. Recent models demonstrate that it 
can be explained by l arge variations in the abundance of stable 
iron group elements ( Wo osley et alJl2007l) with the dominant 
cause for diversity lik ely residing in the explosion mechanism 
dMazzali et al.ll2007l) . 

One additional variable is the metallicity of the WD core, 
which yields excess neutrons relative to protons due to the 
isotope 22 Ne. This is usually expressed as 

i 

where A,- and Z,- are the nucleon number and charge of 
species i with mass fraction X,. The neutronization is crit- 
ical for setting the prod uction of the neutron-rich isotopes 
dThielemann et al.l[l9 86). If no weak interactions occur dur- 
ing the explosion, the mass fraction of 56 Ni produced is 
simply X( 56 Ni) = 587 e - 28, assuming 56 Ni and 58 Ni are the 
only burning products (Timmes et al. 2003). The neutroniza- 
tion also affe cts the explosive burni ng, including the laminar 
flame speed (Chamu fak et al.ll2007al) . However, the metallic- 
ity range of progenitors is not large enough to account for the 
full SNe la diversity (see §4), making it critical to explore all 
factors that determine Y e . 

A potential neutronization site is the convective car- 
bon burning core that is active for ~ 1000 years prior 
to the explosion. The hydrostatic evolution associated 
with this simmering phase terminates when the core tem- 
perature is sufficiently high that burning is dynamical 
dNomoto et al.ll984tlWooslev & Weaver! 198ft [Wooslev et all 
I2004t IWunsch & Wooslevl l2004t iKuhlen et al] 120061). and a 
flame commences (iTimmes & Woosleyll 19921) . We show here 
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that protons from the 12 C( 12 C,/?) 23 Na reaction during simmer- 
ing capture on 12 C, and that subsequent electron captures on 
13 N and 23 Na decrease Y e . This process continues until either 
the explosion occurs or sufficient heavy elements have been 
synthesized that they capture the protons instead. 

In ^2] we present simmering WD core models and summa- 
rize the reaction chains that alter Y e . We find that one proton 
is converted to a neutron for each six 12 C nuclei consumed 
for burning at p < 1.7 x 10 9 g cm" 3 . At densities above this, 
an additional conversion occurs from an electron capture on 
23 Na. Hence, the Y e in the core depends on the amount of 
carbon burned during simmering and the density at which it 
occurs, which we quantify in §3. We find that neutronization 
during simmering dominates for metallicities Z/Z Q < 2/3. 
We conclude in i2]by discussing the observations and noting 
where future work is needed. 

2. NEUTRON PRODUCTION DURING SIMMERING 

Thermally unstable burning begins when the energy gen- 
eratio n rate from carbon fusion, e, exceeds neutrino cool- 
ing {Romoto et al. 1984). The thin solid lines in Figure 
[D show the range o f ignition curves for X( 12 C) = X( 16 0) 
(TYakovlev et all 120061) with the middle line the nominal 
current best. The carbon fuses via 12 C( 12 C,/?) 23 Na and 
12 C( 12 C,a) 20 Ne with branching ratios of 0.44 and 0.56, re- 
spectively. At "early" times the liberated protons capture onto 
12 C, while at "late" times enough heavy elements ( 23 Na or 
23 Ne) have been produced that they capture the protons in- 
stead. 

We treat the evolution during the simmering phase as a 
series of hydrostatic models consisting of an adiabatic con- 
vective core and an isothermal surface at 10 8 K. As long 
as the convection zone is well described as an adiabat this 
is sufficient for resolving the thermal structure without the 
need to explicitly solve the energy transfer equation. These 
assumptions become weaker once the central temperature is 
Tc ^ 7 x 10 8 K, so that burning occurs sufficiently quickly 
that there is considerable ener gy generation within a convec - 
tive eddy overturn timescale (Garc ia-Senz & Woosleyl[l995l) . 
The energy generation does come into play because it sets the 
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heating timescale, % = c p T c /e, where c p is the specific heat 
of the liquid ions (we use linear mixing and the Coulomb en- 
ergy from Chabrier & Potekhin 1998), nearly given by the 
classical Dulong-Petit law c p w IKbI [iiin p , where /i,- is the ion 
mean molecular weight. Since we evaluate tj, using the cen- 
tral conditions it is a lower limit since it should include the 
entire heat capacity of the convective region (Piro & Chang 
2007; see related discussion for neutron stars in Weinberg et 
al. 2006). In this way, for a given thermal profile there is a 
well-defined heating timescale, which connects our stationary 
models to the true time evolution. The thick dashed lines in 
Figure [T]trace out the trajectory of the central temperature, T c , 
and density, p c , for M = 1.35M and M = 1.37Mh (left and 
right, respectively), both using compositions of X( 2 C) = 0.5, 
X( ie O) = 0.48, and X( 22 Ne) = 0.02. These indica te that ^de- 
creases with increasing T c dLesaffre et al. Il2006t |Pk3[2007). 
The thick solid lines show thermal profiles near the end of the 
simmering. 

The simmering phase ends when sub-sonic convection can 
no longer transport the heat outwards because the timescale 
of heating is now less than the convective overturn timescale. 
Since the overturn timescale depends on the integrated energy 
generation rate near the WD center (while we desire a local 
measure of when convection should end for simplicity), we 
assume that this occurs when f/, ~ fd yn = (Gp c )~ 1 ^ 2 , the dy- 
namical timescale. This gives reasonable agreement to other 
more careful calculat ions that find the si mmering phase ends 
when T c w 8 x 10 8 K (IWooslev et al.l2004l) . We plot t h = 10f dyn 
as a dotted line in FigureQ]to indicate where simmering ends, 
since the strong temperature sensitivity of 12 C fusion makes 
this line rather insensitive to the choice of prefactor. If the 
simmering phase ends earlier, it can be considered in the con- 
text of our models by just truncating our results at a slightly 
lower T c . 

2.1. Main Reaction Cycle at Early Times 

At early times, only 12 C, I6 0, or 22 Ne are poten- 
tial proton capture nuc l ei. W e compared these rates us- 
ingJCaughlan & Fowler (1988), including strong screening 
( Salpe ter & van Hornl 119691) . The le O(/?,7) 17 F reaction is 
negligible, whereas resonances in the 22 Ne(p,7) 23 Na reac- 
tion make its rate comparable to I2 C(/?,7) 13 N. However, the 
larger abundance of 12 C means that it captures more pro- 
tons by a factor of (22/12)X( I2 C)/X( 22 Ne) « 40. The fate of 
the synthesized 13 N requires some discussion, as the branch- 
ing amongst the three relevant reactions: 13 N(e~, u e ) 13 C, 
13 N(7,/?) 12 C, 13 N(p,7) 14 0, depends on T, p, and proton mass 
fraction, X p . 

The production of protons is always the rate limiting step, 
so that each proton is almost immediately captured by 12 C. 
This means that we can find X p by balancing the proton pro- 
duction rate from carbon fusion, \n\ 2 ((Tv) 12+12 (where A = 
0.44 is the branching ratio for the reaction 12 C( 12 C,p) 23 Na 
and «i2 is the I2 C number density) with the proton capture 
rate, n p n\2(av) p+ i2, where n p is the proton number density, 

^( 12 C) (<rv) 12+12 .„ 

Ap-A — r . (Z) 

12 (<7V) p+ l2 

This is plotted in the upper panel of Figure [2] The small 
value of X p confirms our equilibrium assumption for the pro- 
ton abundance, and allows us to show that the 13 N(p,7) 14 
reaction is negligible (bottom panel of Fig. |2]i in comparison 
to electron captures. The bottom panel also shows the rates 
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FIG. 1 . — Conditions in the simmering WD. Thin solid lines are the range of 
carbon ignition curves from Yakovlev et al. (2006). The arching thick dashed 
lines show the central trajectory for the simmering core of a 1.35Mq and 
a 1.37Mq WD (left and right, respectively). Thick solid lines are example 
thermal profiles that are composed of an adiabatic convective core connected 
to an isothermal exterior. Each case is taken near the end of the simmering 
phase, which occurs at the dotted line labeled f;, = 10fd yn . Dashed lines are 
crossing points for nuclear reactions described in the text. 
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FIG. 2. — Mass fractions X p and Z(' 3 N) (upper panel) and 13 N destruction 
rates (lower panel), versus temperature. The dotted line at T = 8 X 10 8 K 
divides where I3 N photodisintegrations becomes important. The density is 
3 X 10 9 g cm" 3 and X( 12 C) = 0.5. 

for the reactions 13 N(e", v e ) n C and 13 N(7,p) 12 C, making it 
clear that electron capture dominates for T < 8 x 10 8 K. 

All of the 13 N comes from the protons synthesized by car- 
bon fusion, therefore the 13 N production rate is equal to that 
of protons. Once again, this is the slowest step so that we 
can find the 13 N abundance by balancing n p n\2{av) p +\2 = 
^ ec ( 13 N)«i3, where R ec ( n N) is the 13 N electron capture rate. 
Captures into excited states are unlikely to be dominant, al- 
lowing us to use the measured ft = 4.65 x 10 3 s from the 
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TABLE 1 

Reaction Summary for Main Cycle 



Reaction 


Thermal Energy Release (MeV) 


12 C( 12 C,p) 23 Na 


2.239 


12 C( 12 C,a) 20 Ne 


4.617 


12 C(p,7) 13 N 


1.944 


13 N(<T,!/ e ) 13 C 





13 C(a,n) 16 


2.214 


12 C(n,7) 13 C 


4.947 


23 Na(<T,i/ e ) 23 Ne a 






11 This reaction only occurs when % > ? C c.23 (Fig.[T}- 

ground-state transitions. (Likewise, we use ft = 1.91 x 10 s s 
for 23 Na electron captures in §2.2.) Combining with equation 
(|2|l gives 



13 



X( n N) = X^[X( l2 C)] 



,2 pN a (<tv) 



12+12 



^ec(' 3 N) 



(3) 



where N\ is Avagodro's number, which is shown in the top 
panel of Figure|2] The network is completed by 13 C(a,n) ie O 
and 12 C(n, 7) 13 C, leading to a composition of one each of 13 C, 
ie O, 20 Ne, 23 Na. 

There are two complications. The first is at high T's where 
photodisintegration of 13 N happens faster than the electron 
captures (above the dashed curve labeled by f ec ,i3 = f p h in Fig. 
1). Chemical balance (p+ 12 C^ I3 N+7) is achieved in this 
limit, fixing the proton to I3 N ratio. The 13 N is then slowly 
removed due to electron captures. The electron captures must 
always balance the proton production, so the 13 N abundance 
remains identical to equation (01. Hence, photodisintegration 
adds steps to the reaction chain (and alters the proton density; 
top panel of Fig. 2) but does not modify the conclusion that all 
protons released in 12 C burning lead to 13 N electron capture. 

The second complication is the reaction 23 Na(e", f e ) 23 Ne 
at high densities (p > 1.7 x 10 9 g cm" 3 ). This occurs to the 
right of the dashed line labeled as f/, = f ec 23 in Figure [1] illus- 
trating that these electron captures only take place at certain 
times during the simmering phase, which we account for in 
§3. Electron captures on 13 N would not have time to occur 
above the dashed line labeled f/, = f ec 13 in FigureQ] but this is 
always after the explosion. 

The main reaction cycle is summarized in Table 1. Six 12 C 
are consumed, producing 13 C, I6 0, 20 Ne, and depending on ?/,, 
either 23 Ne or 23 Na. Therefore, during each cycle either one 
or two protons have been converted to neutrons. 

2.2. Late Time Truncation of Neutron Production 

The carbon burning ashes eventually become abundant 
enough to compete with 12 C for proton captures. The rele- 
vant products are 20 Ne, and either 23 Na or 23 Ne. The 20 Ne 
reactions of 20 Ne(p,7) 21 Na and 20 Ne(p,a) l7 F are negligible, 
so we focus on 23 Na and 23 Ne. 

We plot the proton destruction rate for the most relevant 
reactions in Figure [3] as a function of /, the fraction of 12 C 
that has burned. The number density of either 23 Na or 23 Ne is 
taken to be equal to the number density of I2 C burned times 
2/3 x A w 1/3. Circles denote where the rates cross each 
other, which is nearly independent of p. For f/, < ? eCi 23, / must 
exceed w 0.13 before the 23 Na(p, a) 20 Ne reaction becomes 
important, and ends neutronization. 

When th > f e c,23, 23 Ne forms, which has a higher proton 
capture cross section than 23 Na. Burning only / = 0.061 is 
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FIG. 3. — The proton destruction rates for different processes as a function 
of the fraction of 12 C that has burned, /. The number density of 23 Na or 
23 Ne is taken to be /ni2/3. Circles denote crucial places where the proton 
destruction mechanism switches. 



enough that Ne(p,«) Na becomes the primary proton sink. 
Since this makes Na and liberates a neutron, the new 23 Na 
may electron capture again to make 23 Ne, and the reaction 
chain 23 Ne(p,n) Na(e, f f ) 23 Ne can repeatedly occur, making 
many free neutrons. Competing with this process are other re- 
acti on chains that burn these ne utrons (for example, see Table 
6 of lArnett & Thielemann! 19851) . In this regime, it is difficult 
for us to estimate all the key nuclear reactions that will take 
place. Although further neutronization is pos sible, we cannot 
follow this without a full reaction network (Chamulaket al. 
2007b). Such calculations must also be coupled to a realistic 
model for the core temperature evolution (such as what we 
present here). 

3. MAXIMUM NEUTRONIZATION ESTIMATES 

We set j] as the number of protons that are converted to 
neutrons for every six 12 C consumed, so that rj = 2 (77 = 1) for 
th > f ec,23 (th < f e c,23 X where we approximate A « 0.5. The 77 = 
2 case is an upper limit since in parts of the convection zone 
where p < 1 .7 x 10 9 g cm" 3 the 23 Na does not electron capture 
(and 23 Ne that is mixed to lower densities by convection may 
decay). The total neutronization is measured via 



1 2X( 56 Fe) X( 22 Ne) „r)X( n C) 



56 



22 



6x 12 



(4) 



which includes the initial 56 Fe and 22 Ne content. Neutroniza- 
tion halts either when the WD explodes or when freshly made 
heavy elements compete for protons (Fig. 3). 

In the case of competition from fresh heavy elements, trun- 
cation at high densities occurs when / = 0.061 with 77 = 2. The 
maximum change in Y e is therefore 



-8.5 x 10" 



i AT 2 Q 
0.5 



(5) 



A similar limit pertains at lower densities. One way to ex- 
ceed this limit in the high density case is if additional reaction 
chains occur (see §2.2). We show AT c max as a dot-dashed 
line in Figure 4, in comparison to the AF e 's that result from 
X( 22 Ne) = 0.007 and 0.02 (dotted lines). By coincidence, the 
maximum effect of neutronization during simmering is com- 
parable to that associated with a solar metallicity. 

The other possible limiter of neutronization is the onset of 
the explosion. The reactions in Table 1 show that Q « 16 
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MeV is released as thermal energy when six carbon nuclei 
are burned. 2 If we let E c be the total thermal content that is 
within the convective core with respect to the initial isother- 
mal WD, this implies a change AY e = -r]E c m p /QM c in a con- 
vective core of mass M r , 



AK = -6.5 x 10 



Mr. 



-4V f 

2 10 49 ergs M c 



(6) 



For this to compete with the 22 Ne contribution, a total energy 



49 2X( 22 Ne)M c 
1.4 x 10 ergs -- 



t] 0.02 M Q ' 



(7) 



or 7 x 10 15 ergs g must be released prior to the explosion. 

Simmering ends whe n dynamical burning is triggered, re- 
quiring T c 8 x 10 8 K (IWooslev et al.l|2004i) . If the burning 
occurred within a single zone with the specific heat of §2, 
then reaching this T c would require fs 1.3 x 10 16 ergs g , in 
excess of that implied by equation (0. Of course, in reality 
the convective zone extends outward, so that little mass is at 
T c . To accurately determine the resulting neutronization, we 
construct hydrostatic WD models consisting of fully convec- 
tive cores as described at the beginning of ^2] We consider 
isothermal temperatures of either 10 8 K or 2 x 10 8 K. At any 
given mom ent there is a well defined M c dLesaffre et al. 12006c 
lPiro|20 07), and we evaluate the current thermal content by in- 
tegrating the specific heat relative to the initially isothermal 
WD, 



Mr 



c P [T(M)-Tj]dM, 



(8) 



where 7} is the isothermal WD temperature. In this way we 
use our time independent models to find the fraction of carbon 
that must have burned, /, and the associated AY e as T c and M c 
increase with time. We assume no neutrino losses and thus all 
ks 16 MeV of thermal energy contributes to heating. 

In Figure|4]we summarize the results of these calculations. 
In each case, the slope of AY e shows a break at the transi- 
tion from rj = 2 (?/, > ? ec ,23) to 77 = 1. This break occurs later 
for more massive WDs (Fig. [TJ, thus these have more neu- 
tronization during simmering. Increasing the isothermal tem- 
perature decreases M c , so that it takes less burning to reach a 
given T c . These fully integrated models make it clear that sub- 
stantial neutronization occurs prior to the explosion. In com- 
parison to the AY e from 22 Ne, simmering effects dominate if 
X( 22 Ne) < 0.013 orZ/Z Q < 2/3. This thwarts the occurrence 
of high Y e SNe la in low metallicity progenitors. 

4. CONCLUSION AND DISCUSSION 

We have found that significant neutronization of the WD 
core occurs throughout the simmering stage of carbon burning 
until the onset of the explosio n. If substantial energ y is lost to 
the convective Urea process dLesaffre et al. II2005L and refer- 
ences therein), then the neutronization is truncated by proton 
captures onto freshly synthesized heavy elements (resulting 
in eq. @). The main consequence is a uniform "floor" to the 
amount of neutronization that dominates over the metallicity 
dependent contribution for all progenitors with Z/Z Q < 2/3. 

Given the likely significance this has for SNe la, more work 
needs to be done. In particular, at high ignition densities, 

2 Not all of this energy will always go into heating up the core. For exam- 
ple, if the convective Urea process is operating, then it will take more energy 
(and more carbon burning) to get to sufficient temperatures. 
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FIG. 4.— Thick dashed (M = 1.35M Q )and solid (M= 1.37Af ) lines show 
the change of Y e as a function of the current central temperature, T c , due 
to carbon burning during simmering. In each case we consider isothermal 
temperatures of 10 8 K and 2 X 10 s K (as labeled). The dot-dashed line is the 
maximum change possible if proton captures on carbon are only limited by 
competition from freshly made heavy elements (eq. (5)). Dotted lines bracket 
the change in Y e expected for X( 22 Ne) = 0.007-0.02. 



heavy element electron captures and a full reaction network 
are needed to follow the resulting diverse collection of ele- 
ments (see the discussion in §2.2). The convective Urea pro- 
cess is another complication we have not addressed. In prin- 
ciple, if more energy is lost to neutrinos then more burning 
(and thus more neutronization) is required to make the burn- 
ing dynamical. A ssessing this will necess itate coupling a full 
nuclear network (Cham ulak et al.l l2007bl) to convective cal- 
culations. Such models would accurately determine r\ rather 
than simply setting it to 1 or 2. 

In closing, we highlight some important features exhib- 
ited by recent observations of SNe la. It is clear that the 
amount of 56 Ni produced in SNe la has a dynamic range 
(0.1- IMq) larger than can be explained by metallicity or sim- 
mering neutronization. However, since an intriguing trend is 
the prevalence of 56 Ni rich events in star-forming regions it is 
interesting to quantitatively explore how large the observed 
discrepancy is. Using the SNLS sample of Sullivan et al. 
(2006), Howell et al. (2007) found that the average stretch is 
s = 0.95 in passive galaxies (e.g. E/SO's) and s = 1.05 in star- 
forming galaxies. Using Jha et al's (2006) AMis(B)-s rela- 
tion and Mazzali et al.'s (2007) relation between AMis(B) and 
56 Ni mass we get 0.58M Q (s = 0.95) and 0.72M Q (s = 1.05). 
Hence, amongst the large diversity, there is a tendency for 
SNe in star-forming galaxies to produce w 0.1 3M© more 56 Ni 
than those in large ellipticals. 

Since the SN la rate scales with mass in ellipticals and 
star formation rate in spirals (Mannucci et al. 2005; Scan- 
napieco & Bildsten 2005; Sullivan et al. 2006), SNe from 
passive galaxies in the SNLS survey are from more massive 
galaxies than the SNe in star-forming galaxies (Sullivan et al. 
2006). Using the mass-metallicity relation of Tremonti et al. 
(2004), our integration of the separate samples in Sullivan et 
al. (2006) yield average 12 + log(0/H) = 8.87 in active galax- 
ies and 9.1 in ellipticals (solar value is 12+log(0/H) = 8.69). 
Due to the increase of N/O at high metallicities (Liang et 
al. 2006), the SNe in ellipticals have twice as much 22 Ne 
content as those in spirals. From the relation of Timmes et 



NEUTRONIZATION DURING SN la SIMMERING 



5 



al. (2003), this implies sa 5% less 56 Ni, whereas the ob- 
served decrement is > 15%. Explaining the observed decre- 
ment would require a contrast of AX( 22 Ne) ps 0.06, or nearly 
3 times solar. Although we have found that simmering en- 
hances neutronization, the effect is not great enough (AY e max 
would give the same change in neutronization as doubling a 
solar metallicity), and a diverse set of core conditions would 
still be required. A large enhancement could be present in 
the core if substantial gravitational separation h ad occurred 
dBildsten & HalfeOOUlDelove & Bildstenll2002l) . yet convec- 
tive mixing during simmering will reduce it based on the frac- 



tion of the star that is convective. For a convection zone that 
extends out to M©, the resulting 22 Ne enhancement would be 
at most ps 30%. 
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